        r01 = k01(i) * (s(i,HI)) * (s(i,de))
        r02 = k02(i) * (s(i,HII)) * (s(i,de))
        r61 = 0*k01(i) * (s(i,DI)) * (s(i,de))
        r62 = 0*k02(i) * (s(i,DII)) * (s(i,de))
        r03 = k03(i) * (s(i,HeI) / 4.0) * (s(i,de))
        r04 = k04(i) * (s(i,HeII) / 4.0) * (s(i,de))
        r05 = k05(i) * (s(i,HeII) / 4.0) * (s(i,de))
        r06 = k06(i) * (s(i,HeIII) / 4.0) * (s(i,de))
        r07 = k07(i) * (s(i,HI)) * (s(i,de))
        r08 = k08(i) * (s(i,HM)) * (s(i,HI))
        r09 = k09(i) * (s(i,HI)) * (s(i,HII))
        r10 = k10(i) * (s(i,H2II) / 2.0) * (s(i,HI))
        r11 = k11(i) * (s(i,H2I) / 2.0) * (s(i,HII))
        r12 = k12(i) * (s(i,H2I) / 2.0) * (s(i,de))
        r13 = k13(i) * (s(i,H2I) / 2.0) * (s(i,HI))
        r14 = k14(i) * (s(i,HM)) * (s(i,de))
        r15 = k15(i) * (s(i,HM)) * (s(i,HI))
        r16 = k16(i) * (s(i,HM)) * (s(i,HII))
        r17 = k17(i) * (s(i,HM)) * (s(i,HII))
        r18 = k18(i) * (s(i,H2II) / 2.0) * (s(i,de))
        r19 = k19(i) * (s(i,H2II) / 2.0) * (s(i,HM))
        r21 = k21(i) * (s(i,HI))*(s(i,HI))*(s(i,H2I) / 2.0)
        r22 = k22(i) * (s(i,HI))*(s(i,HI))*(s(i,HI))
        r23 = k23(i) * (s(i,H2I) / 2.0)*(s(i,H2I) / 2.0)
        r50 = 0*k50(i) * (s(i,HII) * s(i,DI)/2.0)
        r51 = 0*k51(i) * (s(i,HI) * s(i,DII)/2.0)
        r52 = 0*k52(i) * (s(i,H2I)/2.0 * s(i,DII)/2.0)
        r53 = 0*k53(i) * (s(i,HDI)/3.0 * s(i,HII))
        r54 = 0*k54(i) * (s(i,H2I)/2.0 * s(i,DI))
        r55 = 0*k55(i) * (s(i,HDI)/3.0 * s(i,HI))
